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^^ ' Abstract. In this paper we first calculate the post-Newtonian gravitoelectric secular rate of the mean anomaly of 

^ ' a test particle freely orbiting a spherically symmetric central mass. Then, we propose a novel approach to suitably 

combine the presently available planetary ranging data to Mercury, Venus and Mars in order to determine, simul- 
taneously and independently of each other, the Sun's quadrupole mass moment J2Q and the secular advances of 
the perihelion and the mean anomaly. This would also allow to obtain the PPN parameters 7 and /3 independently. 
We propose to analyze the time series of three linear combinations of the observational residuals of the rates of 
^L ' the nodes Q, the longitudes of perihelia n and mean anomalies M of Mercury, Venus and Mars suitably built up 

f^ , in order to absorb the secular precessions induced by the solar oblateness and the post-Newtonian gravitoelectric 

forces. The values of the three investigated parameters can be obtained by fitting the expected linear trends with 
straight lines, determining their slopes in arcseconds per century and suitably normalizing them. According to 
the present-day EPM2000 and DE405 ephemerides accuracy, the obtainable precision would be of the order of 
^ ' 10~''-10~® for the PPN parameters and, more interestingly, of 10~® for J20. It must be pointed out that the future 

BepiColombo mission should improve the knowledge of the Mercury's orbit perhaps by one order of magnitude. 
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1. Introduction 

In this paper we will deal with certain Newtonian and post-Newtonian secular effects affecting the nodes fi, the 
longitudes of the perihelia n and the mean anomalies A4 of the Solar System's planets. 

Historically, one of the first classical tests of the Einstein General Theory of Relativity (GTR) was the successful 
explanation of the anomalous secular perihelion advance of Mercury in the gravitational field of Sun (Einstein 1915). 
As we will see, this feature of the planetary motion is strictly connected with the problem of the quadrupole mass 
moment J20 of Sun. In (Pireaux and Rozelot 2003) a theoretical range J20 = (2 ± 0.4) x 10~^ is admitted for the 
Sun's oblateness. For this topic and the interplay between both effects see the recent review (Pireaux & Rozelot 2003) 
and the references therein and (Ciufolini & Wheeler 1995). Basically, the point is the following. In regard to the 
secular orbital motions, to which we are interested here, the Sun's mass quadrupole moment induces classical effects 
which have, qualitatively, the same temporal signature of the relativistic ones. This means that they could corrupt 
the recovery of the genuine post-Newtonian features of motion because they could not be removed from the time 
series without removing the post-Newtonian signal of interest as well. All depends on the precision with which J20 
is known: should the mismodcUing in it induce classical residual precessions larger than the relativistic ones, there 
would be no hope to get a reliable test of relativistic gravity. As we will see later, the same problems could come 
from the iV— body secular precessions. In this paper we propose to disentangle such effects by measuring them in an 
independent way. More precisely, we propose to extend a certain approach used in Earth artificial satellite motion 
analysis to the interplanetary arena in order to single out just certain post-Newtonian (and Newtonian) orbital motion 
features independently of the Parameterized-Post-Newtonian (PPN) framework (Will 1993) which is usually employed 
in testing competing metric theories of gravity. Instead, the current approach consists of testing the post-Newtonian 
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Table 1. Orbital parameters of Mercury, Venus and Mars (http://nssdc.gsfc.nasa.gov/planetary/factsheet/). 
For the Astronomical Unit (A.U.) we use the value 1 A.U.=149597870691 m of the DE405 (Standish 1998) and 
EPM2000 ephemerides (Pitjeva 2001a; 2001b). The angle e refers to the inclination to the ecliptic. 

Planet a [A.U.) e (°) e 

Mercury 0.38709893 7.00487 0.20563069 

Venus 0.72333199 3.39471 0.00677323 

Mars 1.52366231 1.85061 0.09341233 



equations of motion as a whole in terms of the PPN parameters which are determined from multi-parameters fits 
together with other astrodynamical quantities. 

In Tabled the orbital parameters of the inner planets of the Solar System are reported. 

2. The post-Newtonian effects 

2.1. The gravitoelectric effects 

In the framework of the standard PPN formalism the post-Newtonian gravitoelectric acceleration induced by the 
Schwarzschild-like part of the spacetime metric and experienced by a test body freely falling around a static, spherically 
symmetric central mass M is (Soffel 1989) 



GM 



GM 

2(7 + /3) ^{v-v) 



r + 2{l+-i){r-v)v\, (1) 



where G is the Newtonian gravitational constant, c is the speed of light in vacuum, r and v are the position and velocity 
vectors, respectively, of the test body, 7 and (3 are the standard Eddington-Robertson-Schiff PPN parameters (7 — 
/3 = 1 in GTR; in this case eq.l^ reduces to the expression of (Mashhoon et al. 2001)). Note that r here is the standard 
isotropic radial coordinate, not to be confused with the Schwarzschild radial coordinate r = r[l + GM/(2c^r)]^. In 
the usual orbital data reductions it is just r which is employed. 

By considering it as a small perturbation to the Newtonian monopole acceleration it is possible to work out its 
effect on the orbital motion of a test body with the standard perturbative techniques. The secular rate of the argument 
of pericentre is given by the well known formula 

ZnGM 
"'^^ - cMl-e^)''^^' ^^^ 

where 

2 + 27-/3 .„. 

J^GE = 5 , (3) 



a and e are the semimajor axis and the eccentricity, respectively, of the test particle's orbit and n = ^GM ja? is the 
(unperturbed) Keplerian mean motion. The orbital period of an unperturbed, two-body Keplerian ellipse is P = 2t: /n. 
Note that eq.lO is an exact result valid to all order in e. 

As we will show, also the mean anomaly M is affected by a post-Newtonian gravitoelectric secular rate. The Gauss 
perturbative equation for M. is given by 

where J7 and i are the longitude of the ascending node and the inclination of the test particle's orbit to the equator 
of the central mass, respectively, and Afj is the radial component of the perturbing acceleration. In order to obtain 
the post-Newtonian gravitoelectric secular rate of the mean anomaly one must consider the radial component of 
eq-Oi insert it in the second term of the right-hand-side of eq.Q), evaluate it on the unperturbed Keplerian ellipse, 
characterized by r = a(l — e^)/(l -I- e cos/) and by the radial and along-track components of the velocity vector which 
are vji = nae sin f /\/l — e? , vt = na{l + ecos f)/\/l — e^, respectively (/ is the orbiter's true anomaly), multiplying 
it by 

M ^ {l-e'f'^df 

P 27r(l-hecos/)2' ^ '' 
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and integrating over one orbital revolution, i.e. from to 27r. It turns out that the post-Newtonian gravitoelectric 
secular rate of the mean anomaly is, to order 0{e^) 



nGM 
Mg^ 5 — 



We can define 



(2 + 47 + 3/5) (l + y)+(2 + 7)e2 



(6) 



2 + 47 + 3/3 
AiGE = ^ . (7) 

It turns out that the second term of eq.© induces for Mercury an additional shift of 1.789 arcseconds/century (" 
cy^^ in the following), while for the other planets it is of the order of 10^^-10^"* " cy^^. Thus, for them the secular 
rate of the mean anomaly can be written as 

Mg^ ; — /- ^ MGE- (8) 

c ay 1 — e 

However, as we will see later, the mean anomaly of Mercury will not be used in the combined residuals strategy outlined 
in the following. For Mercury, the post-Newtonian gravitoelectric effect induced by the Sun on uj and AA , according to 
GTR, amounts to 42.980 " cy~^ and -127.949 " cy~^, respectively. The most accurate estimate of the gravitoelectric 
perihelion advance seems to be that obtained for Mercury by E.M. Standish in 2000 with the DE405 ephemerides 
(Standish 1998) and reported in (Pireaux and Rozelot 2003). He averaged the Mercury's perihelion evolution over 
two centuries by using the DE405 ephemerides with and without the post-Newtonian accelerations. Standish included 
in the force models also the solar oblateness with J20 = 2 x 10~^, so that the residuals for Mercury accounted for 
the post-Newtonian effects only; the determined shift was 42.98 ± 0.0023 " cy^^. The same approach for the mean 
anomaly (Standish E M private communication) has yielded to —130.003 + 0.0027 " cy~^ (See also Table|31later). Note 
that the quoted uncertainty do not come from direct observational errors. They depend on the fact that in the force 
models used in the numerical propagation many astrodynamical parameters occur (masses of planets, asteroids, etc.); 
their numerical values come from multiparameter fits of real data and, consequently, are affected by observational 
errors. Such numerical tests say nothing about if GTR is correct or not; they just give an idea of what would be the 
obtainable accuracy set up by our knowledge of the Solar System arena if the Einstein theory of gravitation would be 
true. 

In regard to the possibility of constructing time series of planetary mean anomalies, it must be noted that a certain 
effort would be required. Indeed, so far this orbital elements has never been utilized, so that the partials, e.g., should 
be computed (Pitjeva, private communication 2004). 

2.2. The gravitomagnetic Lense-Thirring effect 

Another post-Newtonian secular precession which affects not only ui but also 51 is the Lense-Thirring effect (Lense & 
Thirring 1918) induced by the proper angular momentum J of the central body 

—6GJcosi , , 

2GJ 

^LT = ^2„3(l_e2)3/2'^LT, (10) 

where 

1 + T 

MLT = -^. (11) 

By assuming for the Sun Jq — 1.9 x 10^* g cm^ s^^ (Pijpers 2003), the gravitomagnetic effect on, e.g., Mercury's 
perihelion is of the order of 10"'^ " cy"^ (De Sitter 1916). Such small value is at the edge of the present sensitivity 
(0.002 " cy~^) in determining the Mercury perihelion shift from the ephemerides (cfr. the results obtained by Standish). 
Moreover, it should be considered that, even if future improvements of the obtainable experimental sensitivity in inter- 
planetary ranging allowed to consider the possibility of measuring the Lense-Thirring perihelion advance of Mercury, 
the impact of the systematic errors due to the uncertainties in the solar oblateness would severely limit the realistic 
accuracy obtainable in such demanding measurement. Indeed, by assuming an uncertainty of cj^g = 0.4 x 10^^, the 
error in the secular precession induced by the Sun's quadrupole moment would amount to 245% of the Lense-Thrring 
shift of TT^'^'^'^ (see below for the definition of tt). Another important source of systematic error in the measurement of 
such a tiny effect would be represented by the iV— body classical secular precessions which are of the order of lO'^ — 10'^ 
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" cy~^ (see on the WEB http://ssd.jpl.nasa.gov/elem4Dlanets.html). Indeed, their residual mismodelled part 
could severely bias the recovery of the Lense-Thirring effect. For a recent approach to these problems and possible 
strategies to overcome them, see (lorio 2004). It must be noted that, if, on the other hand, we look at a Lense-Thirring 
test as a way to measure the Sun's angular monientuni by assuming the validity of GTR, a measurement of the solar 
gravitomagnetic field would have a significance only if the obtainable accuracy was better than 10%; indeed, among 
other things, the present-day uncertainty in the Sun's angular momentum J©, which could be measured from the 
Lense-Thirring precessions, is just of the order of 10% (Xu & Ni 1997) in various solar models (Paterno et al. 1996; 
Elsworth et al. 1995) or even less in the framework of asteroseismogyrometry (Pijpers 2003). 

At present, the only performed attempts to explicitly extract the Lense-Thirring signature from the data of orbiting 
masses in the Solar System are due to Ciufolini and coworkers who analyzed the laser-ranged data of the orbits of 
the existing LAGEOS and LAGEOS II Earth artificial satellites (CiufoHni et al. 1998). A 20-30% precision level in 
measuring the terrestrial gravitomagnetic field is claimed, but other scientists judge these evaluations too optimistic 
and propose different error budget (Ries et al. 2003). In April 2004 the GP-B spacecraft has been launched. It will 
carry out a very complex and challenging mission which should be able to measure a gravitomagnetic precession of 
the spins of four superconducting gyroscopes (Schiff 1960) carried on board at a claimed accuracy of 1% or better 
(Everitt et al. 2001). 

3. The solar oblateness 

The solar quadrupole mass moment J20 is an important astrophysical parameter whose precise knowledge could 
yield many information about the inner structure and dynamics of our star. A reliable evaluation of J20 still faces 
some controversy: on one side, the theoretical values strongly depend on the solar model used, whereas accurate 
measurements are very difficult to obtain from observations. For all this matter see the recent review (Pireaux & 
Rozelot 2003) and (Rozelot et al. 2004). From an observational point of view, J20 is not directly accessible. In this 
context a dynamical determination of J20, analyzing, e.g., the orbits of the inner planets of the Solar System, is 
interesting because it might be compared with those derived from solar model dependent values of the oblateness. 
However, it is not simple to reach this goal because of the interplay between the effects of the solar quadrupole 
moment with those induced by the post-Newtonian gravitoelectromagnetic forces. Instead of the trajectory of planets, 
it would be possible to infer J20 from accurate tracking of some drag-free spacecraft orbiting within a few radii of the 
solar center. This will be the approach followed by, e.g. the BepiColombo mission (see Section IT!^ . Alternatively, the 
Sun's quadrupole mass moment can be inferred from in-orbit measurement of solar properties, like the SOiiO-MDI 
space-based observations (Armstrong & Kuhn 1999), or from Earth-based observations like those realized, e.g., with 
the scanning heliometer of the Pic du Midi Observatory (Rozelot et al. 2004). 

3.1. The classical precessions induced by the solar oblateness 

For a given planet of the Solar System orbiting the Sun, apart from the classical effects induced by the precession 
of the equinoxes and by the other planets and major asteroids which are routinely accounted for in the ephemerides 
computations (Pitjeva 2001a; 2001b), the oblateness of Sun induces also secular precessions on i7, tt ~ uj + n cosi and 
Ai given by 

3 n. J20 / i?0 



3 nJ20 f^\^ [^ 
'2(l-e2)2 { a J {2'' 



^J2e = -:^n :^{-r] (oSin'«-l), (13) 
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^J^e = I n 2W2 h? (3cos2^-l), (14) 



For Mercury^ the Newtonian precessions due to Sun oblateness, with J20 = 2 x 10 ^, are of the order of 10 ^ " cy ^ 
for il, TT and A4. 

4. The interplay between the solar oblateness and the post- Newtonian precessions 

In Table 121 the relevant parameters of the classical and post-Newtonian secular precessions of the nodes, the perihelia 
and the mean anomalies of Mercury, Venus and Mars are reported. In all the relativistic tests performed up to now by 

^ As pointed out in (Milani et al, 2002), the angle i refers to the inclination between the planet's orbital plane and the fixed 
reference plane of the celestial reference frame; it is not the angle e between the planet's orbital plane and the ecliptic. It turns 
out that i ~ e/2. For Mercury e = 7.00487°. 
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Table 2. Post-Newtonian precessions and coefficients of Newtonian precessions of the node, the perihehon and the 
mean anomaly for Mercury, Venus and Mars in " cy^^. The values ttge: and A^ge are calculated with GTR. The 
coefficients (l,2, t^.2 and AI.2 are d{ilj^Q)/d{J2Q), d{-kj^pj)/d{J2Q) and d{M j^^) / d{J2Q) , respectively. In order to have 
the precessions they must be multiplied by J20 . Note that the result for the mean anomaly of Mercury accounts for 
the correction of order O(e^) which, instead, can be neglected for the other planets. 



Planet JIge ttge A^ge f^.2 i^.-z M.2 

Mercury 42.981 -127.949 -126878.626 126404.437 123703.132 

Venus 8.624 -25.874 -13068.273 13056.803 13056.504 

Mars 1.351 -4.035 -980.609 980.353 976.067 



analyzing the perihelia advances only of the inner planets of the Solar System with the radar ranging technique (Shapiro 
et al. 1972; 1976; Shapiro 1990) it has been impossible to disentangle the genuine post-Newtonian gravitoelectric 
contribution of eq.© from the Newtonian precession of ea. (|13() . Indeed, the observational residuals of tt for a single 
planet, built up by suitably switching off the post-Newtonian 0{c^^) terms and the oblateness of Sun in the force 
models of the equations of motion in the orbital processors softwares, account entirely for the post-Newtonian and the 
solar oblateness^ effects. This is a unsatisfactory situation, both if we are interested in testing post-Newtonian gravity 
and if we want to obtain a dynamical, model- independent measurement of J20. Indeed, it is, of course, impossible to 
constraint both the effects if only one perihelion rate is examined one at a time: in recovering one of the two effects we 
are forced to consider the other one as if it was known. Since the post-Newtonian gravitoelectric effect is three orders 
of magnitude larger than that induced by solar oblateness, a determination of the latter by assuming the validity of 
GTR would be affected by a non negligible systematic error induced by the precision to which the post-Newtonian 
pericentre advance is known from other (more or less indirect and more or less biased by other aliasing effects) tests 
(Lunar Laser Ranging, binary pulsars periastron advance^). The inverse situation is more favorable: indeed, if we are 
interested in the post-Newtonian gravitoelectric effect the relative systematic error induced on its measurement by the 
precession due to the solar oblateness amounts to 5 x 10~^ even by assuming for the latter effect a 100% uncertainty'*. 

4.1. The present-day approach to test post-Newtonian gravity 

At this point it may be interesting to clarify what is the current approach in testing post-Newtonian gravity from 
planetary data analysis followed by, e.g., the Jet Propulsion Laboratory (JPL). In the interplay between the real data 
and the equations of motions, which include also the post-Newtonian accelerations expressed in terms of the various 
PPN parameters, a set of astrodynamical parameters, among which there are also 7 and /3, are simultaneously and 
straightforwardly fitted and adjusted and a correlation matrix is also released. This means that the post-Newtonian 
equations of motion are globally tested as a whole in terms of, among other parameters, 7 and /?; no attention is paid 
to this or that particular feature of the post-Newtonian accelerations. The point is that the standard PPN formalism 
refers to the alternative theories of gravitation which are metric, i.e. based on a symmetric spacetime metric. But 
it is not proven that an alternative theory of gravitation must necessarily be a metric one. Moreover, even in the 
framework of the metric alternative theories, the PPN formalism based on 10 parameters is not sufficient to describe 
every conceivable metric theory of gravitation at the post-Newtonian order; it only describes those theories with a 
particularly simple post-Newtonian limit. One would, in principle, need an infinite set of new parameters to add to 
the standard ten parameter PPN formalism in order to describe the post-Newtonian approximation of any a priori 
conceivable metric theory of gravity (Ciufolini 1991; Ciufolini & Wheeler 1995). 



^ Also if the solar oblateness is included in the force models, the related uncertainty induces a corresponding systematic error 
in the recovered post-Newtonian effect. Fortunately, it is small; by assuming o-j^^ = 0.4 x 10~ , it amounts to 0.01%. 

^ Note that the binary pulsars periastron measurement should not be considered as a test of relativistic gravity because the 
masses of the binary system are not known (Stairs et al. 1998); they can be obtained by assuming the validity of GTR. The 
Lunar Laser Ranging (LLR) measurements do not allow to single out uniquely the gravitoelctric pericentre advance from the 
other post-Newtonian features of motion of the Earth-Moon system (Nordtvedt 2001). Indeed, LLR tests the post-Newtonian 
equations of motion as a whole. Recently, it has been proposed to measure the relativistic gravitoelectric perigee advance of the 
terrestrial LAGEOS II satellite (lorio et al. 2002; Lucchesi 2003), but, up to now, the test has not yet been performed. 

* Note that the impossibility of disentangle the gravitoelectric and Lense-Thirring effects would not seriously affect the recovery 
of the gravitoelctric precession: indeed, the bias induced by the gravitomagnetic effect on the gravitoelectric shift amounts to 
0.004% only for Mercury. 
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4.2. The possibilities opened by the future missions 

Concerning the possibility of disentangle the effects of the solar oblateness from those of the post-Newtonian gravi- 
toelectric force, it is stated that the future space mission BepiColombo^ of the European Space Agency (ESA) will 
provide us, among other things, with a dynamical, model-independent and relativity-independent measurement of J20 
by measuring with high precision the nodal motion of Mercury (Milani et al. 2002; Pireaux & Rozelot 2003) which 
is not affected by the post-Newtonian gravitoelectric force. The claimed accuracy would amount to ctjjq = 2 x 10^^ 
(Milani et al. 2002). However, such evaluation refers to the formal, statistical obtainable uncertainty only. Indeed, 
the residuals of the Mercury's node would account, to a certain level of accuracy, for the Lense-Thirring precession 
as well. By considering such effect as totally unmodelled in the force models, its impact on the measurement of J20 
would induce a 8 x 10~^ systematic error. The formal, statistical accuracy for 7 and /3 is evaluated to be of the order 
of 2 X 10~^ (Milani et al. 2002). Also the ESA astrometric mission GAIA^ should measure, among other things, the 
solar quadrupole mass moment by analyzing the longitudes of the ascending nodes of many minor bodies of the Solar 
System. The obtainable accuracy for 7 is of the order of 10~^-10~^ (Vecchiato et al. 2003). The ASTROD mission 
should be able to measure J20 with a claimed accuracy of the order of 10^* or, perhaps, 10^^ — 10~^° (Ni et al. 2004). 
The claimed obtainable accuracy for the PPN parameters is 4.6 x 10^^ for 7 and 4 x 10~^ for /3. Further improvements 
may push these limits down to 10~^— 10~^. The recently proposed LATOR mission should be able to measure, among 
other things, 7 to a 10~® accuracy level and J20 to a 10~^ level (Turyshev et al. 2004). The solar orbit relativity test 
SORT (Melliti et al. 2002), which would combine a time-delay experiment with a light deflection test, should allow to 
reach a 10~^ accuracy in measuring the PPN parameters 



5. The experimental accuracy in planetary radar ranging 

Concerning the present-day accuracy of the planetary radar ranging, the radar itself is accurate well below the 100 
m level (Standish 2002). The problem, however, comes from the fact that the surfaces of the planets have large 
topographical variations. They are modeled in different ways. For Mercury, spherical harmonics and some closure 
analysis (comparing values when two different measurements reflect off from the same spot on the surface) have been 
done (Anderson et al. 1996) in DE405. For Venus, a topographical model, which comes from (Pettcngill et al. 1980), 
has been used. For Mars, closure points can be used. Closure points are pairs of days during which the observed points 
on the surface of Mars are nearly identical with respect to their longitudes and latitudes on Mars. Since the same 
topographical features are observed during each of the two days, the uncertainty introduced by the topography may 
be eliminated by subtracting the residuals of one day from the corresponding ones of the other day. The remaining 
difference is then due to only the ephemeris drift between the two days. The closure points for Mars have a priori 
uncertainties of about 100 m or less when the points are within 0.2 degrees of each other on the martian surface. Of 
course, for Mars, there is also the spacecraft ranging - far more accurate than the radar: Viking Landers (1976-82), 10 
m; Mars Global Surveyor and Odyssey (1999-2003), 2-3 m. However, correction for Mars topography is possible not 
only by using closure points (in this method no all observations may be used) , but with help of modern hypsometric 
maps and by the representations of the global topography with an expansion of spherical functions (Pitjeva 2001b). It 
must now be noted that our knowledge of the orbital motion of Mercury should improve thanks to the future hermean 
missions Messenger (see on the WEB http://messenger.jhuapl.edu/ and http://discovery.nasa.gov/messenger.html), 
which has been launched in the summer 2004 and whose encounter with Mercury is scheduled for 2011, and, especially^, 
BepiColombo. A complete error analysis for the range and range-rate measurements can be found in (less & Boscagli 
2001). According to them, a full 5-way link to the main orbitcr will be adopted. A multi-channel combination of the 
data will allow to remove most of the measurement errors introduced, in a single channel, by the plasma. As a result, 
a two orders of magnitude improvement in the Earth-Mercury range should be possible, also avoiding the problems 
related to the surface topography. According to a more conservative evaluation by E.M Standish (Standish, private 
communication 2004), improvements in the Mercury's orbital parameters might amount to one order of magnitude, 
i.e. the tens of meters level. 
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See on the WEB http://astro.estec.esa.nl/BepiColombo/. Present ESA plans are for a launch in 2010-2012. 

See on the WEB http://astro.estec.esa.nl/GAIA/. Present ESA plans are for a launch in niid-2010. 

While the spacecraft trajectory will be determined from the range-rate data, the planet's orbit will be retrieved from the 
range data (Milani et al. 2002). In particular, the determination of the planetary centre of mass is important to this goal which 
can be better reached by a not too elliptical spacecraft's orbit. The relatively moderate ellipticity of the planned 400 x 1500 km 
polar orbit of BepiColombo, contrary to the much more elliptical path of Messenger, is, then, well adequate. 
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Table 3. Present-day accuracy in determining the node, perihelion and mean anomaly secular rates of Mercury, Venus 
and Mars according to DE405 (Standish 1998) and EPM2000 (Pitjeva 2001a; 2001b) ephemerides. The figures, in " 
cy~^, represent the formal, statistical errors. Realistic errors should be 10 times larger, at least. While the results by 
Standish come from the mathematical propagation of the nodes, the perihelia and the mean anomalies (Standish, E. 
M. private communication) evolution with and without post-Newtonian terms (with j — (3 — 1) and their average 
over a time span of two centuries, the results by Pitjeva are based on real data. The figure for Mercury has been 
obtained in 2001, while the other ones have been determined subsequently (Pitjeva, personal communication 2004). 
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Mars 



0.000182 
0.000006 
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0.0023 
0.0414 
0.0014 



0.0027 
0.0414 
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0.0086 
0.1037 
0.0001 



6. The proposed approach for disentangle the solar oblateness and the post-Newtonian effects 

Let us write down the following equations® 
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where jn^,^!"''*, 571^^^"°* and SAd^lf^'^^ are the observational residuals^ of the rates of the nodes, the longitudes of 
the perihelia and the mean anomalies of Mercury, Venus and Mars. It is intended that all kind of data (optical and 
radio) would be used. The residuals should be built up by purposely switching off the solar quadrupole moment and 
the post-Newtonian gravitoelectric accelerations in the force models (or leaving in them some default values to be 
subsequently adjusted according to the present strategy) of the orbital processors. Then, the so obtained observational 
residuals would entirely (or partly, if some default values are left in the force models) adsorb just the investigated 
secular effects and other post-Newtonian short-periodic features, i.e. not averaged over one orbital revolution of the 
planet under consideration. In respect to the latter point, it should be noted that the residuals of the mean anomalies 
would account, e.g., also for the indirect effects on the mean motions n through the perturbations in the semimajor 
axes a 



Ar 



3 GM 



Aa. 



(18) 



There are no secular perturbations induced on a by the other planets. If the classical short-periodic effects on a would 
be of relatively no importance because they would be included in the force models at the best of their accuracy, this is 

* Note that in the right-hand-sides of egs. 1161 - 1171 1 also the mismodelled parts of the classical A^— body precessions should 
have been included. Since they are small with respect to the gravitoelectric effects of interest, as we will show later, we can 
neglect them in the calculations for setting up our combiantions. This, of course, does not mean that the experimental residuals 
do not account also for them. 

^ Here we speak about residuals of Keplerian orbital elements in a, strictly speaking, improper sense. The Keplerian orbital 
elements are not directly observable: they can only be computed. The basic observable quantities are ranges, range-rates and 
angles. Here we mean the differences between the time series of the node got from a given observed orbital arc and the time 
series of the node got from a propagated orbital arc with a given force, which we are interested in, switched off in the force 
models. The two time series share the same initial conditions. 
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not the case for the post-Newtonian ones. The gravitoelectric field induces no secular variations on a, as the classical 
planetary perturbations. The short-term shift on a can be calculated from eq.QJ and the Gauss equation for the 
perturbed rate of semimajor axis 

-^= n = [^flesin/-H^T(l + ecos/)], (19) 

at n^/\ — e^ 

where At is the along-track component of the perturbing acceleration. It amounts to 

AflGE = ^7^^^^[14(cos/o - cos/) + 10e(cos2 / - cos^ /q)] + 0{e^). (20) 

c^(l — e'^Y 

For Venus their nominal amplitudes are of the order of 100 m; for Mars they amount to 2 km; ca. H18|) would yield 
periodic variations whose nominal amplitudes would be of the order of 0.2-1 " cy~^ for Venus and Mars, respectively. 
However, such harmonic signatures could be fitted and removed from the data over sufficiently long time spans. Indeed, 
since we are interested in the gravitoelectric secular trends on tt and A^ it should be possible, in principle, to construct 
the residuals by using orbital arcs longer than the sidereal revolution periods of the planets to be used. Then, all the 
high-frequency perturbations would not affect their time series which should, instead, be characterized by the secular 
parts of those Newtonian and post-Newtonian features present in the real data but (partly) absent in the force models 
of the equations of motion in the orbital processors. However, in regard to the possibility of constructing accurate 
time series of observational residuals many years long the following observations must be kept in mind (Standish 
2002). The planetary motions are perturbed by the presence of many asteroids whose masses are quite poorly known. 
Furthermore, it is not possible to solve for the asteroid masses, other than for the biggest few, because there are too 
many of them for the data to support such an effort. As a result, the ephemerides of the inner planets, especially Mars, 
will deteriorate over time; the ephemerides have uncertainties at the 1-2 km level over the span of the observations 
and growing at the rate of a few km/decade outside that span. On the other hand, it must also be noted that the 
sidereal orbital periods of Venus and Mars amount to 224.701 and 686.980 days, respectively. 

We can consider eas. ((T5|l -(fT7 |l as three systems of algebraic linear equations in the three unknowns^° >/20, Mge and 
z/GE- Their solutions can be written as 

SMlT + ^'iSM^^r = MGE (M^r + c[Mlr') , (21) 



+ c:Sn^^T = ^ge(^^ 



A-:rMars _l _ r 'More ,, / ■ Mars i_ „ 'Merc 



where 



■rcf-)Mars oMcrc o Venus oVcnusf-jN 



^--l oVcnusoMara _oA4ars o Venus i . i O A^ i ^ 02 oVcnusoAlars ^fSMarsoVcnus (.llO^U 

"lT "N-body "lT "N-body "lT "N-body "lT "N-body 

^ -TW'f '^'VTW'^™"" = -0.07475, (^^) 

= -^Mars/^Morc ^ _0.00775, 

and 

j^Merc _^ ^^fjvenus _^ ^2j=2^ais ^ -32808.8816 " cy-\ 

■A^ge"" + c'i^Ige"'" = -2.1007 " cy-\ (23) 

Ti-Mars ^ c'l TT^f ^^ = 1.0176 " Cy"! 

The first equation of ea. (|^ comes from ea. l(T3|) solved for J20; it allows to obtain J20 independently of the post- 
Newtonian Lense-Thirring and classical iV— body secular precessions which would represent the major sources of 
systematic errors. The second equation of eq.(j2J comes from ea. lfTB|) solved for /xge; it cancels out the secular 
precessions due to J20. The third equation in eq.|(2U comes from ea. ((T7|l solved for z/qe; it cancels out the secular 
precessions due to Jiq. In regard to the impact of the residual mismodelled classical TV— body precessions on the 

This approach is analogous to that employed in the LAGEOS-LAGEOS 11 Lense-Thirring experiment in the gravitational 
field of Earth (Ciufolini 1996; Ciufolini et al. 1998) 
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second and the third combinations of eq. H21|) , they should not induce a systematic error larger than the observational 
one (see below) because the expected values of /iqe and i^qe are of the order of unity^^, contrary to J20 which should 
be of the order of 10~^. The adimensional parameters ^2©; Mge and z^ge are estimated by fitting the time series of 
the left-hand-sides of ea. H21|) with straight lines, measuring their slopes, in " cy~^, and, then, by dividing them by the 
the quantities of ea. (|23|l which have the dimensions of " cy~^. Note that the solar quadrupole mass moment would 
not be affected by the indirect effects on n because only the nodes would be used in its determination. Finally, from 
the so obtained values of /ige and z/qe, which are 1 in GTR and in Newtonian mechanics, it is possible to measure 
7 and (3 independently of the solar oblateness and also of each other as 

7 = ^(/^ge + !^ge) - |, 

(24) 

/5 = 5MGE — 5^GE + 5- 

According to the results of Table El it is possible to yield an estimate of the (formal) uncertainty in J20, /ige and vq^ 
as 



0'./2£ 



'I'GE 



4.5 X 10^^ 

6 X 10-^ (25) 

3 X 10"^ 




The evaluation for ct^^e accounts also for the fact that eg. 1)181) and the values of aa of Table VI of (Pitjeva 2001a) 
yield a (formal) uncertainty of 3 x 10~^. Note that the accuracies in i^ge and J20 would be further improved from the 
better knowledge of Mercury's orbit which should be obtained from the BepiColombo mission. The uncertainties in 7 
and (3 would amount to 

(26) 

However, discretion is advised in evaluating the reliability of these results because they refer to the formal, standard 
statistical errors; realistic errors may be also one order of magnitude larger. 

For the most recent determinations of 7 and /3 we have that, according to the frequency shift of radio photons to 
and from the Cassini spacecraft, cr^ = 2.3 x 10~^ (Bertotti et al. 2003). This result, combined with ct^ = 4.5 x 10^* for 
the Strong Equivalence Principle (SEP) violating parameter 77 = 4/3 — 7 — 3 (Nordtvedt 1968a; 1968b) from the most 
recent analysis of LLR data, yields u^ = 1 x 10"'* (Turyshev et al. 2003). Conservative results from previous LLR 
analyses yield realistic uncertainties^^ (Miiller et al. 1998) cr-y = 5 x 10"^ and cr^ = 4 x 10"'^. Recent radar ranging 
measurements to planets yield (Pitjeva 2003) cr-y = 1 x 10^'* and ap = 1 x lO^"'. However, as pointed out in (Pitjeva 
2003), these uncertainties are formal standard deviations; realistic error bounds may be an order of magnitude larger. 

7. Conclusions 

In this paper we have explicitly worked out the post-Newtonian gravitoelectric secular rate of the mean anomaly of 
a test particle freely orbiting a spherically symmetric central object. Moreover, we have outlined a possible strategy 
for determining simultaneously and independently of each other, the solar quadrupole mass moment J20 and two 
parameters z^ge and /iGE which account for the post-Newtonian gravitoelectric secular shifts of the perihelion and 
the mean anomaly, respectively, of planets. They are in Newtonian mechanics and 1 in the General Theory of 
Relativity. They could be expressed in terms of the standard PPN 7 and (3 parameters which could, then, be determined 
independently of each other as well. The usual approach employed in the ephemerides data reductions tests post- 
Newtonian gravity theories as a whole straightforwardly in terms of the PPN parameters involving the simultaneous 
fit of many more or less correlated astrodynamical parameters among which there are also 7, /3 and J20. In this case, 



'^^ Let us quantitatively discuss this point. By using the results for the observed centennial rates Q and ro released at 
http://ssd.jpl.nasa.gov/elem_planets.html in, say, the left-hand-side of the third combination of eg. 1211 1 it is possible 
to obtain for it a nominal A'^— body shift of 539.6036 " cy^^. The uncertainty in the A'^— body precessions lies mainly in the Gm 
of the perturbing planets, among which Jupiter plays the major role. Now, the relative uncertainty in Jupiter's Gm is of the 
order of 10"^ m"* s~^ (Jacobson 2003); then, a reasonable estimate of the order of magnitude of the mismodelled part of the 
A'^— body shift should be 1 x 10~® " cy~^. This figure must be divided by 1.0176 " cy~^ yielding a relative error in uqe of the 
order of 10"^ Cfr. with eg.ll^. 

^^ The SEP parameter uncertainty in (Miiller et al. 1998) is cr^ = 9 x 10""*. 
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instead, we propose to analyze the time series of three suitably linear combinations of the residuals of r2, tt and M of 
Mercury, Venus and Mars built up in order to single out just certain selected Newtonian and post-Newtonian orbital 
effects. This can be done by setting purposely equal to zero (or to some default values to be subsequently adjusted 
with the proposed strategy) the orbital effects of interest in the force models of the equations of motion. By suitably 
choosing the length of the orbital arcs it would be possible to account for the secular terms only. The coefficients of 
the combinations would make each of such combinations sensitive just to one orbital effect at a time, independently of 
the other ones. By fitting the experimental residual signals with straight lines, measuring their slopes in " cy^^ and 
suitably normalizing them would yield the values of J20, Mge and vge- The obtainable formal accuracy would be of 
the order of 10~^ — 10~^ for 7 and P and, more interestingly, 10"^ for J20. These estimates should be improved by 
the future more accurate tracking data of Mercury's orbit from the BepiColombo mission. 
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